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The aim of this brief review is to present a case study of how astrophysics data can be used to 
get bounds on Lorentz- violating parameters. For this purpose, a particularly simple Lorentz- 
violating modification of the Maxwell theory of photons is considered, which maintains gauge 
i-G ' invariance, CPT, and renormalization. With a standard spin-| Dirac particle minimally 

Mh| coupled to this nonstandard photon, the resulting modified-quantum-electrodynamics model 

ai involves nineteen dimensionless "deformation parameters." Ten of these parameters lead to 

Q_) ' birefringence and are already tightly constrained by astrophysics. New bounds on the remain- 

ed , ing nine nonbirefringent parameters have been obtained from the inferred absence of vacuum 

Cherenkov radiation in ultrahigh-energy-cosmic-ray (UHECR) events. The resulting astro- 
physics bounds improve considerably upon current laboratory bounds and the implications of 
this "null experiment" may be profound, both for elementary particle physics and cosmology. 
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1 Introduction 

One of the fundamental questions of physics is the following: does space remain smooth as one 
qq | probes smaller and smaller distances? 

A conservative limit on the typical length scale £ of any nontrivial small-scale structure of 

space comes from particle-collider experiments which see no need to change Minkowski space- 
k>( , time: 

LEP/Tevatron: £ < 1(T 18 m « he/ (200 GeV) . (1) 

Yet, astrophysics (or, more specifically, "astroparticle physics") provides us with very much 
higher energies. Here, we intend to present a case study of what astrophysics can do, provided 
the underlying physics is well understood. References will be primarily to research papers. 

The proposed case study proceeds in three steps. First, we discuss the phenomenology of 
a simple photon-propagation model. Second, we obtain bounds on the model parameters from 
ultrahigh-energy cosmic rays (UHECRs). Third, we consider the theoretical implications of our 
results. The contribution concludes with a brief outlook. 

2 Phenomenology 

2.1 Model 

Consider the following action for a Lorentz-violating (LV) deformation of quantum electrody- 
namics (QED): 

S'modQED = -SmodM + SstandD > (2) 



with a modified-Maxwell term [1, 2], 
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modM 



d 4 x ( - \ (»W + k^p°) {d^A u (x) - d u A^x)) (d p A a (x) - d a A p (x))^j , (3a) 



and the standard Dirac term for a spin-Tj particle with electric charge e and mass M, 
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standD 



d 4 x ^(x)(^(id fl - eA p (x)) - M)ijj(x) 



(3b) 



Here and in the following, natural units are used with c = h = 1, but, occasionally, c or h are 
displayed in order to clarify the physical dimension of a particular expression. The fundamental 
constant c now corresponds to the maximum attainable velocity of the Dirac particle or, more 
importantly, to the causal velocity from the underlying Minkowski spacetime with Cartesian 
coordinates (x M ) = (x°,x) = (ct,^ 1 , x 2 ,x 3 ) and metric g pv {x) = 7] pu = diag(+l, —1, —1, —1) . 

The real dimensionless numbers n^^" 7 in (3a) are considered to be fixed once and for all, 
which makes the model Lorentz noninvariant. However, the action (2) is still gauge-invariant, 
CPT-even, and power-counting renormalizable [3]. Clearly, these properties make the model 
worthwhile to study theoretically and to verify/falsify experimentally. 

In the modified-Maxwell term (3a), K pvpa is a constant background tensor with the same 
symmetries as the Riemann curvature tensor and a double trace condition n^ u pu = 0, so that 
there are 20 — 1 = 19 components. 

As ten birefringent parameters are already constrained at the 10 -32 level [4] by spectro- 
polarimetric measurements of distant astronomical sources, the model can be restricted to the 
nonbirefringent sector by making the following Ansatz [5]: 
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for a symmetric and traceless matrix H pv with 10 — 1 = 9 components. 

Hence, there are nine LV deformation parameters K pu to investigate. It turns out to be 
useful to rewrite these parameters as follows: 
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with a single independent parameter k 00 for the spatially isotropic part of W and eight inde- 
pendent parameters 8H pi/ for the nonisotropic part. Finally, we can express these parameters in 
terms of the so-called standard-model-extension (SME) parameters [4]: 
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where three parity-odd parameters determine an antisymmetric traceless 3x3 matrix (k _| 
and five parity-even parameters a symmetric traceless 3x3 matrix (/5 c _) mn . 
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Figure 1: Left: classical spacetime-foam manifold. Right: Feynman diagram for vacuum Cherenkov radiation. 

2.2 Possible spacetime origin 

We have already mentioned that the modified-Maxwell model has attractive properties (gauge 
invariance, CPT invariance, and renormalizability) . Still, it would be better if we had a concrete 
example of how the modified-Maxwell model might arise from an underlying theory. 

Precisely this is accomplished by a calculation [6] of the propagation of standard photons and 
standard Dirac particles over a classical spacetime-foam manifold (Fig. 1-left), which reproduces 
a restricted, isotropic version of model (2): 

2K tr = -a 2 F, 6^ = 0, (7) 

in terms of the quadratic coefficient of the modified photon dispersion relation to be given 
shortly. 

In fact, the simplest possible classical spacetime-foam model has identical "defects" (equal 
spatial and temporal size b) embedded with random orientations in Minkowski spacetime (equal 
spatial and temporal average separation 7). For this type of manifold, modified dispersion rela- 
tions for protons (p) and photons (7) have been obtained in the large-wavelength approximation: 

u? v = M 2 cph 2 + c 2 p k 2 + 0(k A ), (8a) 

(1 + a 2 F) c 2 p k 2 + (54 PI 2 ) c 2 p k A + 0(fc 6 ) , (8b) 



^ 



with wave number k = |k| = 2ir/\, effective defect on/off factors 02, 04 G {±1> 0}, effective 
defect size b, and effective defect excluded- volume factor, 

F ee (6/I) 4 , (9) 



which is assumed to be less than 1. Remark that the coefficient of the quadratic term in (8a) is 
c 2 



simply defined as c 2 and that F appears in both the quadratic and quartic terms of (8b) [the 



inclusion of F in the quartic photon term defines the meaning of the squared length b 2 ] . 

Calculations of specific space-time-foam models [6] give the effective parameters (with tilde) 
in terms of the underlying spacetime parameters (with bar) : 

b = f3b, 7= A 7, a 2 = -l, a 4 = l, (10) 

for positive constants (3 and A of order unity. From (9) and (10), the restricted (isotropic) model 
(7) is defined by the fundamental length scales of the classical space-time-foam model considered, 
namely, the typical defect size b and the average separation 7. [Note that this classical space- 
time-foam model, in its simplest form, involves a preferred frame of reference, as the assumed 
equality of spatial and temporal defect sizes makes clear.] In principle, spacetime defects can 
also generate CPT-violating terms in the effective photon action [7, 8], but the corresponding 
physical effects are expected to be suppressed by at least one power of the fine-structure constant 
a = e 2 /4vr « 1/137. 



2.3 Vacuum Cherenkov radiation 

In certain Lorentz-violating photon models, the decay process p — ► p7 is allowed. This process 
is similar to that of Cherenkov radiation in material media [9, 10, 11, 12, 13] but now occurs 
already in vacuo due to the modified photon propagation and may, therefore, be called "vacuum 
Cherenkov radiation" [14, 15, 16, 17]. For the modified-QED model (2), the decay process 
P ~> PI has been studied classically by Altschul [18] and quantum-mechanically (Fig. 1-right) 
by Kaufhold and the present author [19]. 

The radiated-energy rate of a point particle with electric charge Ze, mass M > 0, momentum 
q, and ultrarelativistic energy E ~ c \q\ is given by [19] 
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~ —mE 2 /n, (ii) 



with a direction-dependent coefficient £ > and a threshold energy [18, 19] 

<esh = „ ro „ -7^ ^^7^1 + O (M 2 c 4 ) , (12) 



for LV parameters \n ta/ \ <C 1 and ramp function 

R[x] = (x + \x\)/2. (13) 

Note that (11) becomes infinite in the classical limit h — ► 0, which traces back to the fact that 
quantum mechanics provides a frequency cutoff that renders the total radiated energy finite; see 
Refs. [11, 12, 13, 19] for further discussion and references. 

Exact tree-level results have been obtained [20, 21] for the restricted modified-QED model 
(labeled "modQED, isotropic case") with ao = (4/3) k 00 = 2/5 t r > and 8fi& iv = 0, which is 
precisely the model derived in the previous subsection. Setting again h = c = 1, the radiated- 
energy rate for a spin-^ Dirac point particle (charge Ze, mass M, and energy E = A/|q| 2 + M 2 
above threshold) is given by 
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x h(al + 4a + §)E 2 - (2 + a ) U{l + a )M 2 + 2(3 + 2a ) j|— ^ £ y^ 2 - M 2 J I. (14) 
The high-energy expansion of (14) for fixed parameters ao and M reads 
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, (15) 

which displays the quadratic behavior of (11) for a constant coefficient £ = (7/24) oq + O(qq). 
From (14), one also obtains the exact tree-level threshold energy: 

7-,modQED, isotropic case "**■ P, | 777 fin\ 

^thresh = ~f= V 1 + "0/2 , (16) 



which reproduces (12) for a small positive value of «o — 2%- Incidentally, the numerical value 
of ao cannot be too large, as the factors y/2 — a$ i n (14) make clear. 

As mentioned before, the radiated-energy rates (11) and (14) have been derived for point 
particles. For realistic particles (protons, nuclei, and photons), the partonic content can be ex- 
pected to significantly modify the prefactors of the radiation rates but not the energy thresholds 
(12) and (16), which follow from energy-momentum conservation [model (2) violates Lorentz 
invariance but maintains spacetime translation invariance]. Note that calculations of standard 
classical Cherenkov radiation from different types of charge distributions give different radiation 
rates, depending on the details of the distributions, but unmodified energy thresholds, at least 
for the case of static charge distributions in the rest frame and a nondispersive medium (cf. 
Sec. 7.2 of Ref. [13]). a The bounds on the LV parameters (6) derived in the next section rely 
only on the energy threshold (12) and are therefore reliable. 

3 UHECR bounds 

3.1 Cherenkov threshold condition 

The basic idea [14, 15] consists of three steps: 

• if vacuum Cherenkov radiation has a threshold energy ^thresh (*0, then UHECRs with 
-E-prim > -^thresh cannot travel far, as they rapidly radiate away their energy, 

• this implies that, if an UHECR of energy E pr \ m is detected, its energy must be at or below 
threshold, 

-Eprim < -Ethresh(ft), (17) 

• the last inequality gives, using expression (12), an upper bound on the LV parameters [19], 

R[27i tT + 2 5K ^ prim + 5^ k ql rim q'; Tim ] < (M p 2 rim c 4 )/E 2 piim , (18) 

with primary energy £p r i m , flight direction q pr j m , and rest mass M pr ; m as input. 

Remark that the Cherenkov threshold condition (18) is only effective for certain combinations 
of deformation parameters k^ u and flight directions q, namely, if there is a positive argument 
of the ramp function R defined by (13). Hence, there must be a sufficient number of UHECR 



parameters (6). 



events with primary energies -E , p rim and flight directions q pr j m , in order to bound all deformation 



3.2 Bounds on the LV photon parameters k^ v 

The last analysis of Refs. [22, 23] is based on the following 29 selected UHECR events: 27 events 
from Auger-South [24], 1 event from Fly's Eye [25], and 1 event from AGASA [26]. The last two 
events are taken in order to provide a better coverage of the northern celestial hemisphere (later, 
they can perhaps be replaced by Auger-North events). The relevant data of these 29 events are 
given in Table 1, where the uncertainties in the energies are of the order of 25% and those in 
the pointing directions of the order of 1 deg (see the original references for further discussion) . 
Adopting a conservative value for the unknown primary mass (see below), it is this relatively 



a The proton is a complicated dynamic object (as is the photon) and "vacuum Cherenkov radiation" is not 
quite the same as standard Cherenkov radiation in a material medium. For this reason, we have embarked on 
a parton calculation [21] of the proton radiated-energy rate for the isotropic case (7) with only parameter K tr 
nonzero and positive. Preliminary results indicate that the proton radiated-energy rate is suppressed for energies 
E € [-Ethresh, /-Ethresh], with £Wesh given by the point-particle result (16) and / ~ 1.25 (based on the results of 
Table 1 in Ref. [16]). 



Table 1: Selected UHECR events from Auger (2004-2007) [24], Fly's Eye (1991) [25], and AGASA (1993) [26]. 
Shown are the arrival time (year and Julian day), the primary energy E [in EeV = 10 18 eV], and the arrival 
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large energy uncertainty which dominates the error budget of the bounds to be presented in this 
subsection and the next. 

The 29 primary energies and flight directions from Table 1 (flight directions being the oppo- 
site of arrival directions) give 29 inequalities for the deformation parameters from the Cherenkov 
threshold condition (18), where we set M pT [ m = 56 GeV/c 2 . [A significant fraction of these pri- 
maries may well be protons [24], but we prefer to take a conservative value for M pr i m . The used 
value of 56 GeV/c 2 is even larger than the mass of an iron nucleus 56 Fe.] The 29 inequalities 
then give the following two— a Cherenkov bounds [23] on the nine isolated SME parameters (6) 
of the nonbirefringent modified-Maxwell model (2)— (4): 



(ij) G {(23), (31), (12)} 
(kl) G {(11), (12), (13), (22), (23)} 
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(19a) 
(19b) 
(19c) 



for the Sun-centered celestial equatorial coordinate system in Cartesian coordinates. A single 
UHECR event suffices for bound (19c) and the 148 EeV Auger event from Table 1 has been 
used, which has a reliable energy calibration. 

Based on a 212 EeV Auger event [27] and setting M pr i m = 52 GeV/c 2 , a new result [21] 
gives the following two-sided bound at the two-<r level on a single universal isotropic parameter: 
—2 x 10~ 19 < Kj* r mv < 6 x 10~ 20 , where the lower bound arises solely from partonic effects. Here, 
"universal" means that the same LV parameter K tr applies to all gauge bosons of the Standard 
Model, as might be expected from a spacetime-foam model as discussed in Sec. 2.2. 

It is important to realize that the Cherenkov bounds (19abc) only depend on the measured 
energies and flight directions of the charged cosmic-ray primaries at the top of the Earth at- 
mosphere. As noted in Refs. [15, 16], the travel length from vacuum Cherenkov radiation (if 
operative) would be of the order of meters rather than megaparsecs. Hence, we only need to 
be sure of having observed a charged primary traveling over a distance of a kilometer, say, in 
order to apply the Cherenkov threshold condition (18). [The previous discussion applies, strictly 
speaking, only to event energies above 1.25 -©thresh ( as mentioned in Ftn. a), whereas cosmic rays 
with energies barely above ©thresh would still have to travel over astronomical distances.] 



For comparison, the current laboratory bounds are as follows (with selected references): 

• direct bounds [28, 29] on the three parity-odd nonisotropic parameters in k q+ at the 10~ 12 
level; 

• direct bounds [28, 29] on the five parity-even nonisotropic parameters in /? e _ at the 10~ 14 
to 10" 16 levels; 

• direct bound [30, 31] on the single parity-even isotropic parameter n tr at the 10 -7 level; 

• indirect bound [32] on /?tr at the 10 -8 level from the measured value [33, 34] of the electron 
anomalous magnetic moment a e = (g e — 2)/2 . 

Interestingly, the UHECR Cherenkov bounds (19) are the strongest where the laboratory bounds 
are the weakest. A case in point is the isotropic parameter K tr , which is difficult to constrain by 
laboratory experiments (having no sidereal variations) but easy to constrain by a single UHECR 
event (the flight direction being irrelevant). 

The "leverage factor" for the UHECR Cherenkov bounds is, according to (18), given by 
(Ep V i m /M pT - im c 2 ) 2 and is of the order of 10 18 for M prim ~ 50 GeV/c 2 and -Ep rim ~ 50 EeV. 
With further information on selected UHECR events (e.g., the shower-maximum atmospheric 
depth A max used already in Ref. [22]), it may be possible to reduce M pr \ m by a factor 10 to 
-^prim ~ 5 GeV/c 2 and, with more and more events becoming available, it may be possible to 
select on a typical energy £p r i m ~ 150 EeV, thereby increasing the Cherenkov leverage factor 
by a factor 10 3 to a value of the order of 10 . For completeness, the leverage factor for the 
astrophysics bounds [4] on the birefringent parameters is given by L/X, which is of the order of 
10 32 for a source distance L ~ 10 26 m and typical wavelength A ~ 10 -6 m. Both leverage factors 
make clear what the power of astrophysics can be, provided the physics is well understood. 

3.3 Bounds on the LV photon parameter b 2 

From the 148 EeV Auger event in Table 1, we also get a one-a bound [6] on the general coefficient 
of the quartic photon term in (8b): 

\a A Fl 2 \ < (2x KT 35 m) 2 , (20) 

based on the parton analysis of Ref. [16] but rescaled to M pr i m = 56 GeV/c 2 and £p r i m = 
148 EeV. In fact, the right-hand side of (20) is of order (he 3 M prim / E 2 Thn ) 2 . Note that the above 
bound is two-sided, whereas genuine Cherenkov radiation would only give a one-sided bound (the 
phase velocity of the electromagnetic wave must be less than the maximum attainable velocity 
of the charged particle). The reason for having a two-sided bound (20) is that, in addition to 
Cherenkov radiation, another type of process can occur. The relevant process is proton break-up 
p — ► pe + e~ (that is, pair-production by a virtual gauge boson), which then gives the "other" 
side of the bound [16, 21]. 

It may be of interest to mention that the non-observation of primary photons by the Pierre 
Auger Observatory [27, 35] has also been used to obtain (modulo some assumptions) a tight 
one-sided bound [36] on the quartic term of a modified photon dispersion relation. 

In the analysis of Refs. [16, 36], it is taken for granted that the modified photon dispersion 
relation corresponds to a consistent theory with, for example, microcausality and unitarity (cf. 
the discussion in Refs. [37, 38]). This may very well be the case, provided that the cubic term 
is absent from the modified dispersion relation, leaving a quartic term as the first higher-order 
term to be considered (cf. the discussion in Refs. [6, 39]). But it is also clear that a real 
understanding of potential Lorentz-violating effects can only come from considering a complete 
theory, even if it is only an effective theory such as the one discussed in Sec. 2.2. 



As remarked in an earlier review [40], bound (20) disagrees by many orders of magnitude 
with a possible "quantum-gravity" effect [41] in a gamma-ray flare from Mkn 501 as observed 
by the MAGIC telescope [42]. Most likely, the observed time dispersion has an astrophysical 
origin. Still, this type of measurement suggests how, in principle, astrophysical data could give 
more than just upper-bounds on the possible small-scale structure of space. 

The length scale on the right-hand side of (20) is extraordinarily small, but this inequality 
becomes less dramatic if the dimensionless number F on the left-hand side is also small. Now 
this is precisely what has been found in the model calculation leading up to (8), where the 
quartic photon coefficient has an extra reduction factor F which can be interpreted as the 
defect excluded- volume factor (9) entering the quadratic photon coefficient. Taking a value 
F = 1.5 x 10~ 19 , just consistent with the previous result (19c), bound (20) becomes 

b 2 < (5x 10~ 26 m) 2 , (21) 

which remains small compared to the current laboratory bound (1). 

4 Theoretical implications 

In the previous section, we have established two types of bounds on Lorentz violation in the 
photon sector. First, a combined one-u bound [4, 23] was obtained on the nineteen Lorentz- 
violating deformation parameters of the modified-Maxwell model (2): 

\^ vpu \ <3x 10 -18 , (22) 

where, for the sake of argument, the "one-sided" Cherenkov bound on the isotropic parameter 
Ktr has also been made "two-sided" [as mentioned in the paragraph starting a few lines under 
(19c), there is a new two-sided bound at the 10~ 19 level]. 

Second, restricting to the isotropic model (7)— (10) with a possible spacetime-foam origin, 
one-o" Cherenkov-type bounds [6, 16, 21] have been obtained for combinations of the effective 
defect size b and separation I: 

F = (b/T) 4 < 1.5 x 10~ 19 , (23a) 

b < 5 x 10~ 26 m « he/ (4 x 10 9 GeV) , (23b) 

where the particular parameters b and I are really defined by the modified dispersion relations 
(8ab). For simplicity, we focus our discussion on this last case with a single quadratic photon 
coefficient F and a single quartic photon coefficient Fb 2 , both of which can be interpreted in 
terms of a simple spacetime-foam-like structure. 

Bound (23b) is remarkable compared to what can be achieved with particle accelerators on 
Earth (recall that the proton beam energy of the Large Hadron Collider is 7 x 10 3 GeV). As 
it stands, bound (23b) may be easily satisfied by a quantum-gravity theory with length scale 
fcianck = y/hG/c 3 « 1.6 x 10~ 35 m « he/ (1.2 x 10 19 GeV). But for TeV-gravity models [43, 44] 
it would be hard to understand why the effective length scale b for photon propagation in the 
three-dimensional world would be reduced by a numerical factor of the order of 10 6 compared to 
the nonperturbative gravity scale L grav = hc/E grav ~ hc/(4TeV). In fact, the quartic coefficient 
of the photon dispersion relation would need an additional factor of the order of 10~ 12 or less. 
A priori, there is no reason to expect the quartic photon coefficient to be extraordinarily small; 
see Ref. [6] for a heuristic discussion based on so-called Bethe holes [45]. 

Bound (23a) is even more interesting as it implies that a single-scale (b ~ I) classical space- 
time foam is ruled out altogether. This result holds, in fact, for arbitrarily small values of the 
defect size b , as long as a classical spacetime makes sense. In the context of effective theories 



with Lorentz invariance violated at an ultraviolet scale A, a similar result holds that strong LV 
effects can be expected at low energies [46], which have not been seen experimentally. 

The conclusion is, therefore, that Lorentz invariance remains valid down towards smaller and 
smaller distances, which answers in part the question posed at the beginning of the Introduction. 
This conclusion would hold down to distances at which the classical-quantum transition occurs, 
which may happen for distances of the order of /pianck ~ 10~ 35 m or perhaps for distances given 
by an entirely new fundamental length scale [47]. 

5 Outlook 

Astrophysics data (in particular, results from UHECRs) show that a hypothetical quantum 
spacetime foam must have "crystalized" to a classical spacetime manifold which is remarkably 
smooth, as quantified by the defect excluded-volume factor (b/l) 4 < 1CP 19 <C 1 and Lorentz- 
violating parameters \n iMV( " T \ < 10~ 18 <C 1. The result applies to any theory of quantum space- 
time, be it Matrix-theory [48, 49] or loop quantum gravity [50]. 

The outcome is like having a "null experiment" and there is an analogy with the well- 
known Michelson-Morley experiment [51]: theory foresees physical effects which are not found 
by experiment. 

This suggests the need for radically new concepts, similar in depth to the "relativity of 
simultaneity" introduced by Einstein [52]. For our problem, a first small step may be the 
realization that precisely Lorentz invariance is crucial at the high-energy frontier and that a new 
type of conserved relativistic "charge" can play an important role for the flatness of spacetime 
by resolving the so-called cosmological constant problem [53, 54, 55]. It remains to be seen if 
this is a step in the right direction, but experiment has, at least, provided theory with a base 
camp for the long climb towards the quantum origin of spacetime. 
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